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ABSTRACT 

We study non-linear structure formation in high-resolution simulations of Early Dark Energy 
(EDE) cosmologies and compare their evolution with the standard ACDM model. In Early 
Dark Energy models, the impact on structure formation is expected to be particularly strong 
because of the presence of a non-negligible dark energy component even at very high red- 
shift, unlike in standard models that behave like matter-dominated universes at early times. In 
fact, extensions of the spherical top-hat collapse model predict that the virial overdensity and 
linear threshold density for collapse should be modified in EDE model, yielding significant 
modifications in the expected halo mass function. Here we present numerical simulations that 
directly test these expectations. Interestingly, we find that the Sheth & Tormen formalism for 
estimating the abundance of dark matter halos continues to work very well in its standard 
form for the Early Dark Energy cosmologies, contrary to analytic predictions. The residuals 
are even slightly smaller than for ACDM. We also study the virial relationship between mass 
and dark matter velocity dispersion in different dark energ y cosmologies, findi ng excellent 
agreement with the normalization for ACDM as calibrated bv lEvrard et al.l (l2008l) . The earlier 
growth of structure in EDE models relative to ACDM produces large differences in the mass 
functions at high redshift. This could be measured directly by counting groups as a function 
of the line-of-sight velocity dispersion, skirting the ambiguous problem of assigning a mass 
to the halo. Using dark matter substructures as a proxy for member galaxies, we demonstrate 
that even with 3-5 members sufficiently accurate measurements of the halo velocity dispersion 
function are possible. Finally, we determine the concentration-mass relationship for our EDE 
cosmologies. Consistent with the earlier formation time, the EDE halos show higher concen- 
trations at a given halo mass. We fin d that the magnitu de of the difference in concentration is 
well described by the prescription of lEke et"an (1200 ll) for estimating halo concentrations. 
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1 INTRODUCTION 

Arguably the most surprising result of modem cosmology is that 
all matter (including both atoms and non-baryonic dark matter) ac- 
counts for only a quarter of the total energy density of the Uni- 
verse today, while the rest is contributed by a dark energy field. In 
1999, observati ons of type la supemovae by t he Supernovae Cos- 
mology Project jRiess et al]| 19991 : lRiessll2004l) and the relative ac- 
curate measurements of the distances to this objects jPerlmutterl 
ll999l : lKow alski 2008) demonstrated that the expansion of the Uni- 
verse is accelerated today; there hence exists a mysterious force that 
acts against the pull of gravity. Nowadays, the inference that this is 
caused by dark energy can be made with significant confidence, as 
the observational evidence has further firmed up. In fact, we have 
good reason to believe that we live in a flat universe with an upper 
limit of ilm ^ 0.3 for the matter density today, based on cosmic 
microwave ba ckground measureme nts and a host of other observa- 
tional probes jKomatsu et al.ll2008l . e.g.). These observations yield 



a consistent picture, the so-called concordance cosmology, and are 
in agreement with predictions of the inflationary theory. 

The physical origin of dark energy is however unknown and a 
major puzzle for theoretical physics. A nagging outstanding prob- 
lem is that most quantum field theories predict a huge cosmological 
constant from the energy of the quantum vacuum, up to 120 orders 
of magnitude too large. There is hence no simple natural explana- 
tion for dark energy, and one has to be content with phenomenolog- 
ical models at this point. Two proposed forms of dark energy are the 
cosmological constant, a constant energy density filling space ho- 
mogeneously, and scalar fields such as quintessence. In particular, 
'tracking quintessence' models attempt to alleviate the coincidence 
problem of the cosmological constant model. More exotic models 
where the dark energy couples to matter fields or can cluster itself 
have also been proposed. 

In light of the many theoretical possibilities, the hope is that 
future observational constraints on dark energy will enable progress 
in the understanding of this puzzling phenomenon. This requires 
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the exploitation of the subtle influence of dark energy on structure 
formation, both on linear and non-linear scales. As the expected 
effects are generally small for many of the viable dark energy sce- 
narios, it is crucial to be able to calculate structure formation in 
dark energy cosmologies with sufficient precision to tell the dif- 
ferent models apart, and to be able to correctly interpret observa- 
tional data. For example, in order to use the abundance of clusters 
of galaxies at different epochs to measure the expansion history 
of the universe, one needs to reliably know how the cluster mass 
function evolves with time in different dark energy cosmologies. 
Numerical simulations are the most accurate tool available to ob- 
tain the needed theoretical predictions, and they are also crucial for 
testing the results of more simplified analytic calculations. 

In this study, we carry out such non-linear simulations for a 
particular class of dark energy cosmologies, so-called Early Dark 
Energy (EDE) models where dark energy might constitute an ob- 
servable fraction of the total energy density of our Universe at the 
time of matter radiation equality or even big-bang nucleosynthe- 
sis. While in the cosmological constant scenario, the fraction in 
dark energy is negligible at high redshift, in such models the en- 
ergy fraction is a few per cent during recombination and structure 
formation, which introduces interesting effects due to dark energy 
already at high redshift. In particular, for an equal amplitude of 
clustering today, we expect structures to form earlier in such cos- 
mologies than in ACDM. This could be useful to alleviate the ten- 
sion between a low erg normalization suggested by current observa- 
tional constraints from the CMB on one hand, and the observations 
of relatively early reioruzation and the existence of a population of 
massive halos pre sent already at high r edshift on the other hand. 

Recentlv. Bar telmann et al.l ( |2006|) studied two particular EDE 
models, evaluating the primary quantities relevant for structure for- 
mation, such as the linear growth factor of density perturbation, the 
critical density for spherical collapse and the overdensity at virial- 
ization, and finally the halo mass function. In the two models an- 
alyzed, they found that the effect of EDE on the geometry of the 
Universe is only moderate, for example, distance measures can be 
reduced by 8%. Assuming the same expansion rate today, such 
models are younger compared to ACDM. At early times, the age of 
the universe should differ by approximately 5 — 10%. 

However, when Bartelmann et al. (2006) repeated the calcu- 
lation of the spherical collapse model in the EDE cosmology, a 
few nontrivial modifications appeared. The evolution of a homoge- 
neous, spherical overdensity can be traced utilizing both the virial 
theorem and the energy conservation between t he collapse and the 
turn a r ound time (see also Lace y & Cole 1993i : jWang & Steinhardtl 
Il998l) . iBartelmann elal] j2006h obtained the value of the virial 
overdensity as a function of the collapse redshift, translating the 
effect of the early dark energy in an extra contribution to the po- 
tential energy at early times. They found that the virial overdensity 
should be slightly enlarged by EDE, because a faster expansion 
of the universe means that, by the time a perturbation has turned 
around and collapsed to its final radius, a larger density contrast has 
been produced. However, at the same time they found that the lin- 
early extrapolated density contrast corresponding to the collapsed 
object should be significantly reduced. 

These two results based on analytic expectations have a pro- 
nounced influence on the predicted mass function of dark matter 
halos. In EDE models, the cluster population expected from the 
Press-Schechter or Sheth-Tormen formalism grows considerably 
relative to ACDM, as a result of the lowered value of the critical lin- 
ear density contrast 5c for collapse. This effect can be compensated 
for by lowering the normalization parameter ag in order to obtain 



the same abundance of clusters today. In this case, one would how- 
ever still expect a higher cluster abundance in EDE at high redshift, 
due to the earlier growth of structure in this model. 

An open question is whether the EDE really participates in 
the virialization process in the way assumed in the analytic model- 
ing. Similarly, it is not clear whether the excursion set formalism of 
Sheth & Tormen yields an equally accurate description of the non- 
linear mass function of halos in EDE cosmologies as in ACDM. 
Because accurate theoretical predictions for the halo mass function 
are a critical ingredient for constraining cosmological parameters 
(in particular Q,m and SI a) as well as models of galaxy formation, 
it is important to test these predictions for the EDE cosmology in 
detail with numerical N-body simulations. In particular, we want to 
probe whether the fraction / of matter ending up in objects larger 
than a given mass M at some redshift z can be found by only look- 
ing at the properties of the linearly evolved density field at this 
epoch, using the ordinary ST formalism, or whether there is some 
dependence on r edshift, power spectrum or dark energy parameters, 
as suggested bv iBartelmann et al] ( l2006h . 

A further interest in EDE cosmologies stems from the fact that 
for a given erg, the EDE models predict a substantially slower evo- 
lution of the halo population than in the ACDM model. This could 
explain the higher normalization cosmology expected from clus- 
ter studies relative to analysis of the CMB. The value of erg, for 
a given cosmology, provides also a measure of the expected bias- 
ing parameter that relates the galaxy and the mass distribution. The 
early dark energy cosmologies could hence reduce the current mild 
tension between cluster data and the CMB observations. We note 
that halos in cosmologies with EDE are also expected to be more 
concentrated than in ACDM; because the density of the Universe 
was greater at early times, objects that virialized at high redshift are 
more compact than those that virialized more recently. 

Previous numerical simulations of a quintessence component 
with a changing equation of state (EOS) explored two particular 
potentials: SUGRA and Ratra Peebles (RP), which differ because 
RP has a more smoothly decreasin g w and consequently a very 
different evo l ution in the past. Both Linder & Jenkins! | |2003|) and 
iKlvpin et alj l l2003h analyzed the influence of the dark energy on 
the halo mass function in order to extrapolate the abundance of 
structure at different epochs and to compare it with existent theo- 
retical models. They used different numerical codes: the publicly 
available code GAD GET, in the first proj ect, and the Adaptive Re- 
finement Tree code jKravtsov et aljl997h . in the second. They con- 
cluded that the best way to understand which dark energy Universe 
fit the observations best is to look at the gro wth history of halos 
and the evolution of their properties with time. lDolag et al.l ( |2004) 
focused on the modification of the concentration parameter with 
mass and redshift, for the same cosmologies, based on high reso- 
lution simulations of a sample of massive halos. A limited num- 
ber of numerical studies also considered the possibility of a cou- 
pling of the dark en ergy field with dark matter jMainini et alj2003l: 
iMaccio et alj2004l) . 

In this paper, we carry out several high resolution simulations 
of dark energy cosmologies in order to accurately measure the 
quantitative impact of early dark energy on abundance and struc- 
ture of dark matter halos. To this end, we in particular measure halo 
mass functions and evaluate the agreement/disagreement with dif- 
ferent analytic fitting functions. We also test how well the growth 
of the mass function can be tracked with dynamical measure based 
on the velocity dispersion of dark matter substructures, which can 
serve as a proxy for the directly measurable line-of-sight motion of 
galaxies or line widths in observations, and gets around the usual 
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ambiguities arising from different possible mass definitions for ha- 
los. Finally, we also present measurements of halo concentrations, 
and of the relation between dark matter velocity dispersion and halo 
mass. While finalizing this paper, Francis et al. (2008) submitted a 
preprint which also studies numerical simulations of EDE cosmolo- 
gies. Their work provides a different analysis and is complementary 
to our study, but it reaches similar basic conclusions about the halo 
mass function. 

This paper is organized as follows. After a brief introduction 
to the Early Dark Energy models in Section 2, we present the simu- 
lations and also give details on our numerical methods in Section 3. 
In Section 4, we study the mass function of halos for the different 
cosmologies, and as a function of redshift. Then, in Section 5 we 
investigate the properties of halos by studying the virial relation 
between mass and dark matter velocity dispersion, as well as the 
mass-concentration relationship. In Section 6, we consider the ve- 
locity distribution function and prospects for measuring it in obser- 
vations. Finally, we discuss our results and present our conclusions 
in Section 7. 



2 EARLY DARK ENERGY MODELS 

The influence of dark energy on the evolution of the Universe is 
governed by its equation of state, 



p — wpc 



(1) 



A cosmological constant has wa = — 1 at all redshift, while a dis- 
tinctive feature of the Early Dark Energy (hereafter EDE) models as 
well as of other models such as quintessence is that their equation 
of state parameter, Wdc(z), varies during cosmic history. 

Negative pressure at all times implies that the energy density 
parameter will fall to zero very steeply for increasing redshift. If, 
however, we allow the equation of state parameter to rise above 
zero, we can construct models in which ^ldc{z) has a small posi- 
tive value at all epochs, depending on the cosmological background 
model we adopt. While canonical dark energy models with near 
constant behaviour for w do not predict any substantial dark en- 
ergy effect at 2 > 2, in such EDE models the contribution of dark 
energy to the cosmic density can be of order of a few percent even 
at very high redshift. 

We are here investigating this interesting class of models 
which are characterized by a low but non-vanishing dark energy 
density at early times. Note that while the acceleration of the ex- 
pansion of the Universe is a quite recent phenomenon, the dark 
energy responsible for this process could have an old origin. In 
fact, field theoretical models have been c onstructed that ^eneri- 
cally cause such a dynamical behaviour ( Ratra & Peebles 1988"; 
[Wetterich 1988; Ferreira & Jovce 1998; Liddle & Scherrer 1999). 

IWetterichI ( 12004 ) proposed a useful parameterization of a fam- 
ily of cosmological models with EDE in terms of three parameters: 

• the amount of dark energy today, fidc.o (we assume a flat uni- 
verse, so f^m.O = 1 - fidc,0), 

• the equation-of-state parameter wo today, and 

• an average value fidc.c of the energy density parameter at early 
times (to which it asymptotes for z t-^ oo). 

Figure[T]shows the redshift evolution of the equation-of-state 
parameter in the four different cosmologies that we examine in 
this study. As can be noticed, the EDE models approach the cos- 
mological constant scenario at very low redshift. We can compute 
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Figure 1. Equation of state parameter w shown as a function of redshift for 
the four different cosmological models considered in this work. In the two 
early dark energy models EDEl and EDE2, shown with black solid and red 
dotted lines respectively, the value of w today is close to that of ACDM, but 
the amount of dark energy at early times is non vanishing, as described by 
the parameterization (2). 



the equation-of-state parameter for these early dark energy models 
from the fitting formula: 



wi^z') - 

where 
6 = -■ 
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and y = In (1 + 2) = — In a. The parameter h characterizes the 
time at which an approximately constant equation-of-state changes 
its behaviour. 

In Figure |2l we plot the evolution of the matter and energy 
density parameters up to redshift z = 30. The dark energy pa- 
rameter for EDE models evolves relatively slowly with respect to a 
standard ACDM cosmology. In fact, the critical feature of this pa- 
rameterization is a non-vanishing dark energy contribution during 
recombination and structure formation (see also .Doran et al. 200 ih: 



n. 



= — In a„ 



f2do (a) dlna. (4) 

In acq 

For sufficiently low fidc.c, the EDE models reproduce quite well 
the accelerated cosmic expansion in the present-day Universe and 
they can be fine-tuned to agree both with low-redshift observation s 
and CMB temperature fluctuation results dDoran et alj2005l,l2007h . 



3 NUMERICAL SIMULATIONS 

We performed a series of cosmological N-body simulations for 
two early dark energy models 'EDEl' and 'EDE2', which have 
Wo — —0.93 and too ~ —0.99, respectively, and a dark energy 
density at early times of about 10~* (see Tab.[Tll. For comparison, 
we have also calculated a model 'DECDM' with constant equation 
of state parameter equal to w = —0.6, and a conventional ACDM 
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Figure 2. Evolution of the density parameters Q,m{z) (dashed lines) and 
^de(^) (solid line) for the four cosmological models studied in this work. 
At redshift z = 30, the dark energy contribution is orders of magnitude 
higher for EDE models compared with a ACDM cosmology. 



Figure 3. Hubble expansion rate for the models studied in this work. All 
models are normalized with respect to the reference ACDM case. In the 
models EDEl, EDE2, and in the model with constant w, the expansion 
rate of the universe is higher at early times. This has a strong effect on the 
evolution of the growth factor. 
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Table 1. Parameters of the N-Body simulations. The parameter f2(jp de- 
scribes the amount of dark energy at early times, see equation js). This 
value, together with too, the value of the equation state parameter today, 
and f^tjc 0. the amount of dark energy today, completely describes our EDE 
models. 



reference model. We shall refer with these labels to the different 
models throughout the paper. 

In all our models, the matter density parameter today was cho- 
sen as f2m = 0.25, and we consider a flat universe. The Hubble 
parameter k h — _H'o/(100 kms^^Mpc^^) = 0.7 and we as- 
sume Gaussian density fluctuations with a scale-invariant primor- 
dial power spectrum. The normalization of the linear power spec- 
trum extrapolated to 2 = is erg = 0.8 for all our simulations 
in order to match the observed abundance of galaxy clusters to- 
day, irrespective of the cosmology. We also used the same spectral 
index n = —1 throughout in order to focus our attention on possi- 
ble differences due to the dark energy contribution alone. For these 
choices, the models EDEl and EDE2 are almost degenerate, but 
their proximity serves as a useful test for how well differences in 
the results can be detected even for small variations in the EDE pa- 
rameters. This gives a useful illustration on how well one can hope 
to be able to distinguish them in practice and provides realistic data 
for testing the discriminative power of specific statistics. 

For our largest calculations we used 512'^ particles in boxes 
of volume lOO'' ft-'^Mpc^, resulting in a mass resolution of rup = 
5.17 X 10* h~^MQ and a gravitational softening length of e = 
4.2h,-^kpc, kept fixed in comoving coordinates. All the simula- 
tions were started at redshift zinit = 49, and evolved to the present. 
For the si mulations, we adapted the cosmolog ical code GADGET-3 
(based on lSpringel et aljF200lUSpringelll2005l) and the initial con- 
dition code N-GENIC, in order to allow simulations with a time- 
variable equation of state. These simulations can be used to deter- 
mine the mass function also in the high-mass tail with reasonably 



small cosmic variance error, while at the same time probing down 
to interestingly small mass scales. 

In Figure [3] we plot the expansion function of the EDE mod- 
els relative to the ACDM case. We note that the only modification 
required in the simulation code was to update the expression for 
calculating the Hubble expansion rate, which needs to include the 
quintessence component. This term enters in both the kinematics 
and the dynamics of the cosmological models. 

According to the Friedmann equation within a flat universe we 

have 

1/2 

H{a) = Ho + r2dc,o exp ( -3 / [l + w (a)] din a 



(6) 



The density of dark energy changes with the scale factor as: 



f^dc(2:) = f2dc,o exp 



dlna [1 + ui (a) 



instead of simply being equal to f2de,o, as in the usual scenario. For 
ui = — 1, the behaviour of a cosmological constant is recovered. 

If we interpret the modified expansion rate as being due to 
w{z), as defined in equation we find: 



H"" {z) I Hi = fide.O (1 + Z\ 



3+3iu^(z) 



+ n^,o(i + 2) 



(7) 



where 

""'^'^^ i + bin{i + zy 

and b is given by the Eqn. {S}. 

We can see that effectively the EDE models predict the ob- 
served effect of an acceleration in the expansion rate, and this has 
consequences on the global geometry of the Universe. We note that 
the dark energy term in Eqn. ^ just parametrizes our ignorance 
concerning the physical mechanism leading to an increase in ex- 
pansion rate. However, once the dependence of H on the scale 
factor is fixed, the mathematical problem of calculating structure 
growth is then unambiguously defined. 

The evolution of fide, a affects not only the expansion rate of 
the background but also the formation of structures. The primary 
influence of dark energy on the growth of matter density pertur- 
bations is however indirect and arises through the sensitive depen- 
dence of structure growth on the expansion rate of the universe. In 
Figure [4] we show the linear growth factor D divided by the scale 
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Figure 4. Ratio of the growth factor of hnear density perturbations and 
the scale factor a, as a function of a. The four models are described in 
Table [T] The curves are normalized to unity at early times, i.e. we here 
assume that the starting density contrast is the same in the four cosmologies. 
The models EDEl and EDE2 show a significant difference in the growth 
factor evolution even with small energy density at high redshift: structures 
have to grow earlier to reach the same abundance as the ACDM model 
today. 



factor D/a as a function of time for all our models. All curves are 
normalized so that they start from unity at early times. 

In order to rescale the power spectrum of matter fluctuations 
to the redshift of the initial conditions {z — 49 for all simula- 
tions), we introduced in our initial condition code the calculation 
of the growth factor fo r a general equation-of-state as given by 
iLinder & Jenkins! j2003h : 



w (a) D' 
l + X{a)ir 



3 

+ t: 



X{a) D 



2 1 + X{a) 



= 0, 



(9) 



where X (a) is the ratio of the matter density to the energy density: 
f 1 



X{a) 



n 



m,0 



dG,0 



■ exp 



3 / din a'w (a') 



(10) 



and we allowed for a time-dependent equation of state, u;(a). Here 
we define the growth factor as the ratio D = 5 (a) /5 (a^) of the 
perturbation amplitude at scale factor a relative to the one at a,;, 
and we use the normalization condition D(acq) — a^q. 

We can easily see that for very large redshift we recover the 
matter dominated behaviour in the ACDM case: D{a) oc a. On the 
other hand, as expected, the linear growth in the two EDE models 
falls behind the green curve in Fig. [4] implying that they reach a 
given amplitude at earlier times. In fact, the expansion rate in the 
ACDM cosmology is lower than in EDE models, which governs 
the friction term (a /a) in the growth equation 



5 + 2-S~ 4nGpS - 
a 







(11) 



of the perturbations. 

These formulae can be used to derive a suitable expression 
for the r educed linear overdensi ty Sc for collapse expected in EDE 
models l lBartelmann et al.ll2006l) . which in turn suggests that there 
are significant consequences for the process of non-linear structure 
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Figure 5. Comparison of the non-linear power spectra of the four different 
cosmological models studied here. The three panels give results for redshifts 
z = 0, z = 3 and z = 5, from top to bottom. The j/-axes shows the 
dimensionless power = k^P{k) as a function of k computed from the 
dark matter density field using a grid of 512'^ points. All the simulations 
are normaHzed to erg =0.8 for the linearly extrapolated density field today. 
The dashed lines i ndicat e the expected linear power spectra. The prediction 
from lSmith et al.l |2003|) for the ACDM cosmology is shown by the black 
dot-dashed line. 



formation, an expectation that we will analyse later in detail. In 
sum, structures need to grow earlier in EDE models than in ACDM 
in order to reach the same amplitude at the present time. At an 
equal redshift, the initial conditions must hence be more evolved 
in order to produce comparable results today. The DECDM shows 
a behaviour qualitatively similar to EDEl and EDE2 (blue long- 
dashed line). 
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In all our simulations, we have identified dark matter halos us- 
ing two methods: the friends-of-friends (FOF) algorithm with link- 
ing length b = 0.2, and the spherical overdensity (SO) group finder. 
Candidate groups with a minimum of 32 particles were retained by 
the FOF group finder. In the SO algorithm, we first identify FOF 
groups, and then select the particle with the minimum gravitational 
potential as their centres, around which spheres are grown that en- 
close a fixed prescribed mean density A x pcrit, where pcrit is the 
critical density. Different definitions of virial overdensity are in use 
in the literature, and we consider different values for A where ap- 
propriate. The classical definition of NFW adopts A = 200 in- 
dependent of cosmology, while sometimes also A = 200 Q,m is 
used, corresponding to a fixed overdensity relative to the back- 
ground density. Finally, a value of A ~ 178 fim^^ based on a 
generalization of the spherical top-hat collapse model to low den- 
sity cosmologies can also be used. Note however that this may in 
princ iple depend on the dark energy cosmology (Bartelmann et al. 
12006 !), and is hence slightly ambiguous in these cosmologies. 

We have verified the correctness of our implementation of 
early dark energy in the simulation code by checking that it ac- 
curately reproduces the expected linear growth rate in these non- 
standard cosmologies. Recall that rather than normalizing the den- 
sity perturbations of the initial conditions to the same value at the 
(high) starting redshift, we determine them such that they should 
grow to the same linear amplitude today in all of the models. In 
practice, we fix erg, the linearly extrapolated rms fluctuations in 
top-hat spheres of radius 8 /i^^Mpc to the value 0.8 for the epoch 
z = 0. 

In Figure[5l we show measurements of the power spectrum of 
our different models at three different redshifts. While the models 
differ significantly at high redshift, the four different realizations 
show the same amplitude of the power spectrum, at least on large 
scales, at redshift z = Q (left panel). The fluctuations on small 
scales probably reflect the earlier structure formation time in the 
EDE models and the resulting differences in the non-linear halo 
structures. The good agreement of the power spectrum at the end, 
as well as a detailed comparison of the growth rate of the largest 
modes in the box with the linear theory expectation (not shown), 
demonstrate explicitly that the EDE models are simulated accu- 
rately by the code, as intended. 

Note also that the power spectrum measurements show that 
due to the slower evolution of the linear growth factor in EDE 
models, the degeneracy between the models is lifted towards high-z 
since this corresponds to more time for the different growth dynam- 
ics to take effect. Consequently, we expect a different evolution of 
structures back in time. Our main focus in the following will be to 
study the impact of a different equation of state for the dark energy 
upon the mass function of dark matter halos and its evolution with 
redshift. 



4 THE MASS FUNCTION 

In this section we measure the halo abundance at different redshifts 
and compare with analytic fitting functions proposed in the litera- 
ture. Our primary goal is to see to which extent dark energy models 
can still be described by these fitting formulae, and whether there 
is any numerical evidence that suppo rts the higher halo abund ance 
predicted for the EDE cosmologies jBartelmann et alJbOOfil) . We 
will mostly focus on halo mass functions determined with the FOF 
algorithm with a linking length of 0.2, but we shall also consider 
SO mass functions later on. 



In Figure [T] we show our measured halo mass functions in 
terms of the multiplicity function, which we define as 



M dn(M, 

po dluCT" 



(12) 



where po is the background density, n(M, z) is the abundance 
of halos with mass less than A/ at redshift z, and a is the mass 
variance of the power spectrum filtered with a top-hat mass scale 
equal to M. We give results for the cosmological models ACDM, 
DECDM, EDEI and EDE2, plotted as symbols, while the solid 
lines show various theoretical predictions. 

Note that we plot the mass function only in a limited mass 
range in order to avoid being dominated by counting statistics or 
resolution effects. To this end we only consider halos above a min- 
imum size of 200 particles. At the high mass end, individual ob- 
jects are resolved well, but the finite volume of the box limits the 
number of massive rare halos we can detect. We therefore plot the 
mass function only up to the point where the Poisson error reaches 
~ 14% (corresponding to minimum number of ~ 50 objects per 
bin). 

As is well known, the Press & Schechter mass function 
jPress & Schechted [T974h , while qualitatively c orrect, disagrees 
in de t ail with the results of N-body simulation s I Efstathiou et al.l 
Il988l : IWhite et all 1 1 9931 : iLacev & Col j 1 1 994 lEke et al.l 1199^^ 
specifically, the PS formula overestimates the abundance of halos 
near the characteristic mass Af* and underestimates the abundance 
in the high-mass tail. We therefore omit it in our comparison. The 
discrepancy is largely resolved by replacing the spherical collapse 
model of the standard Press & Schechter theory with the refined el- 
lipsoidal collapse model (Sheth & Tormen 1999; Sheth et al. 200 H; 
ISheth & Tormenll2002i) . Indeed, in the top left panel of Figure[6]we 
can see quite good agreement of the Sheth & Tormen mass function 
(ST) with our simulations at 2: = 0. We stress that here the standard 
value of 5c ~ 1.689 for the linear collapse threshold has been used 
irrespective of the cosmological model. Two other well-known fit- 
ting f ormulae are that from Jenkins (central panel, [j enkins et aU 
I2OOIL 'J') and that from Warren (right panel. Iwarren et aU.200a 
'W'), which differ only very slightly in the low-mass range. We 
compare our measurements with these models in the panels of the 
middle and right columns. As we can see from the comparison be- 
tween the solid lines and the numerical data points, the differences 
between the different theoretical models (which only rely on the 
linearly evolved power spectrum at each epoch) and the simulation 
results is very small. 

Figure |7] shows the redshift evolution of the mass function, 
in the form of separate comparison panels at redshifts z — 1 and 
z — 3. While ^t z — the different cosmologies agree rather well 
with each other, as expected based on the identical linear power 
spectra, at redshift z = 1 we begin to see differences between the 
models, and finally at 2 = 3, we can observe a significantly higher 
number density of groups and clusters in the non-standard dark en- 
ergy models. Notice that the model with constant w (blue line) be- 
haves qualitatively rather similar to the EDE models. In each of the 
panels, we include a separate plot of the residuals with respect to 
the analytic fitting functions. This shows that at z = 3 the agree- 
ment is clearly best for the ST formula. 

The differences between the models are most evident in the 
exponential tail of the mass function where it begins to fall off quite 
steeply, in agreement with what is expected from the power spec- 
trum analysis. We can see that, at high-z, replacing the cosmolog- 
ical constant by an early dark energy scenario has a strong impact 
on the history of structure formation. In particular, non-linear struc- 
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Figure 6. Friends-of-friends multiplicity mass functions at z = for the four dark energy models studied here. The solid hnes in each panel represent the 
multiplicity function computed analytically either from the Sheth & Tormen formula (left panel), the Jenkins formula (central panel) or the Warren model 
(right panel). The symbols are the numerical simulation results for ACDM (green), DECDM (blue), EDEl (orange) and EDE2 (red). We consider only halos 
with more than 200 particles and we apply an upper mass cut-off where the Poisson error reaches 14%. In the lower plot of each figure we show the residuals 
between analytically expected and numerically determined mass functions for all models. The differences are typically below 10%. The error bars show 
Poisson uncertainties due to counting statistics for all models. At 2 = 0, the simulation results for all cosmologies are basically identical, which reflects the 
fact that we normalized the models such that they have the same Hnear power spectra today, with a normalization of erg = 0.8. 



tures form substantially earlier in such a model, such that a differ- 
ence in abundance of a factor of ~ 2 is reached already by 2 = 3. 
This underlines the promise high redshift cluster surveys hold for 
distinguishing different cosmological models, and in particular for 
constraining the dynamical evolution of dark energy. 

We now want to assess in a more quantitative fashion the dif- 
ferences between our numerical halo mass functions and the ana- 
lytic fitting functions. In particular we are interested in the question 
whether we can objectively determine a preference for one of the 
analytic models, and whether there is any evidence that the ordi- 
nary mass function formalism does work worse for the generalized 
dark energy models than for ACDM. The latter would indicate that 
the critical linear overdensity threshold 5c needs to be revised for 
EDE models, as suggest ed by the analytic spherical collapse theory 
teartelmann et alj20od) . 

To this end we directly measure the goodness of the fit, which 
we define for the purposes of this analysis as: 




E 



(MF. - MFth.O 



(13) 



where MFth.i are the theoretical values, MFi are the simulations 
results, and we took into account a simple Poisson error in the def- 
inition of the goodness of fit. In Figure [S] we plot this value ex- 
pressed in percent for all simulations when compared with the the- 
oretical formulae of ST (solid line), Jenkins (dotted line) and War- 
ren (dashed line). We cannot identify a clearly superior behaviour 
of any the three fitting functions, at least at this level of resolu- 
tion; the models lie in a strip between approximately 5 and 15% 
error between 2 = and 2: = 5. There is some evidence that the 
ST model does a bit better than the other fitting formulae for the 
ACDM cosmology at high redshift, but the opposite is true for the 
two EDE cosmologies and the Jenkins and Warren functions. 

Interestingly, the overall agreement between simulation re- 
sults and fitting functions is actually slightly worse for ACDM than 
for the non-standard dark energy cosmologies. There is hence no 
tangible evidence that a revision of the mass function formalism 
is required to accurately describe EDE cosmologies. Our finding 
of a universal /(cr) is quantitatively different from the expecta- 
tion based on the analysis of the EDE models bv lBartelmann et al] 



( l2006h . We find that only the different linear growth rate has to 
be taken into account for describing the mass function in the early 
dark energy cosmologies with the ST formalism, but there is no 
need to modify the linear critical overdensity value. To make this 
point more explicit, we show in Figure [TO] the mass function for 
the EDE models and compare it to standard ST (solid lines), and 
to the expectations obtained taking into account a different den- 
sity contrast for EDE models (dashed lines). The predictions in the 
second case are based on the analytic study of iBartelmann et al] 
(12006) and the critical overdensity is proportional to (a)'^^d'=.=f 
(see Eqn. |4j0 Clearly, the proposed modification of Sc actually 
worsens the agreement, both for the halos selected according the 
EOF algorithm (top panel) or defined with respect to the virial over- 
density (bottom panel). 

In the plots we discussed above, we always employed the EOF 
halo finder with standard linking length of & = 0.2 to find the ha- 
los, and the masses were simply the EOF group masses, which ef- 
fectively correspond to the mass within an isodensity surface of 
constant overdensity relative to the background density. As the an- 
alytic mass function formulae have been calibrated with FOE halo 
mass functions, we expect that they work best if the mass is defined 
in this way. However, we may alternatively also employ a different 
mass definition based on the spherical overdensity (SO) approach, 
which allows one to take into account the time-dependent virial 
overdensity A predicted by generalizations of the spherical col- 
lapse model for dark energy cosmologies. In the bottom panel of 
Figure [To] we can see that an even more marked disagreement re- 
sults when we take into account this arguably more consistent halo 
definition. 

To stress this conclusion, in Figure [g] we show the residuals 
of our SO halo mass functions compared with the Sheth & Tormen 
prediction, as a function of redshift and for our different cosmolog- 
ical models, using the same procedure already applied to the EOF 
halo finder results. In this case, the halos were defined as virialized 
regions that are overdense by a variable density threshold equal to 



^ In order to obtain the new values for the critical overdensity it is neces- 
sary to compute the virial overdensity by solving the equation of the gener- 
alized spherical collapse model. 
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Figure 7. Friends-of-friends multiplicity mass functions for the four dark energy models studied here. The evolution towai'ds high redshift is shown in terms 
of results at z = 1 (left column) and at z = 3 (right column). The sohd lines in each plot represent the multiplicity function computed analytically from the 
Sheth & Tormen formula (top row), the Jenkins formula (middle row) and the Warren formula (bottom row). The points are the numerical simulation results 
for ACDM model (green), DECDM (blue), EDEl (orange) and EDE2 (red). We consider only halos with more than 200 particles and we apply an upper mass 
cut-off where the Poisson error reaches 14%. In the lower plot of each figure we show the residuals between analytically expected and numerically determined 
mass functions for all models. The differences are typically below 15%. The error bars show Poisson uncertainties due to counting statistics for all models. 
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Figure 8. Redshift dependence of our goodness of fit parameter (see 
Eqn. ll3t , expressed as percent, computed by comparing the theoretical ex- 
pectation for the multiplicity mass function with the simulation results. All 
cosmological models are compared. The deviations are computed with re- 
spect to the Sheth & Tormen model (solid Unes), the Jenkins et al. (dotted 
lines) and the Warren (dashed lines). 

Ac = IStt^ + 82s - 393;^ (14) 

where x — Qrn{z) — 1, see lBrvan & Normal ( Il998h . This is the 
predicted dependence of A for ACDM, which we used for sim- 
phcity also for the other dark energy cosmologies. As expected, 
we see that the error increases relative to the FOF mass functions, 
with discrepancies of order 10% at 2 = 0. However, there is again 
no evidence that the non-standard dark energy cosmologies are de- 
scribed worse by the ST formalism than ACDM. Also, there is no 
improvement in the accuracy of the fit when we introduce the mod- 
ified linear density contrast for the EDE models. On the contrary, as 
seen by the dotted lines, which represent the theoretical mass func- 
tion (based on Sheth & Tormen) mo dified according to the sph erical 
top hat collapse theory proposed bv lBartelmann et al] ( l2006h . 

Our results thus suggest that the mass function depends pri- 
marily on the linear power spectrum and is only weakly, if at all, 
dependent on the details of the expansion history. This disagrees 
with the expectations from the generalization of the top hat col- 
lapse theory, which are not confirmed by our numerical data. In 
fact, our simulations show that a description of the mass function 
based on the generalized TH calculation is incorrect at the accu- 
racy level reached here. While the dynamic range of our results 
could be improved by increasing the resolution and box-size of our 
simulations, it appears unlikely that this could affect our basic con- 
clusions. Nevertheless, better resolution would be required if one 
seeks to still further reduce the present residuals of order 5-15% 
between the fitting functions of ST, Jenkins or Warren. 



5 HALO PROPERTIES 

5.1 The virial scaling relation 

lEvrard et alj dlOOSh have shown that the dark matter velocity dis- 
persion of halos provides for accurate mass estimates once the re- 
lationship between mass and velocity dispersion is accurately cal- 
ibrated with the help of numerical simulations. They have demon- 



Figure 9. Redshift dependence of our goodness of fit parameter (see 
Eqn. 113) . expressed as percent, computed by comparing the theoretical ex- 
pectation for the multiplicity mass function with the simulation results. All 
models are considered. Here we use the top hat halo mass definition to com- 
pute the mass function from the simulations, and the deviations are com- 
puted with respect to the standard Sheth & Tormen model (dotted-dashed 
lines), and the Sheth & Tormen formula computed from a generalization of 
the top-hat collapse theory (dashed lines). 
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Figure 10. Multiplicity mass function at 2 = 1.5 for the two EDE mod- 
els studied here. We want to highlight that the introduction of a mod- 
ified overdensity motiv ated by the generalized spherical collapse theoiy 
jBartelmann et al.ll2006l) reduces the agreement between the simulation re- 
sults and the theoretical ST mass function. The measured points are in bet- 
ter agreement with the solid line (standard ST model) than with the dashed 
lines (ST modified model), the latter systematically overestimate the halo 
abundance. This holds true both if we consider the halos obtained from 
the FOF halo finder (top panel) or the one obtained by taking into account 
the theoretically motivated virial overdensity with the appropriate SO mass 
definition (bottom panel). 
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Figure 11. The virial scaling relation at the present epoch for primary halos with mass larger than 10 for the four models considered (from left to right 

and top to bottom: ACDM, DECDM, EDEl, EDE2). The red solid line in each plot represents the lEvrard et alj j20Q8) relation, while the blue dashed line is 
our best fit. The triangles are the simulation res ults: we employ a fi xed critical threshold of A = 200 to identify the dark matter halos. The insets show the 
distributions of deviations in In (Tdm around the lEvrard et alj j2008l) fit. 



strated that there exists a quite tight power-law relation between the 
mass of a halo and its one-dimensional velocity dispersion ctdm, 
where 

^ iVp 3 

= - Vj\^ , (15) 

^ i=l j = l 

with Vi,j being the jth component of the physical velocity of par- 
ticle i in the halo, Np is total number of halo particles within a 
radius that encloses a mean overdensity of A = 200 with respect 
to the critical density, and v is the mean halo velocity. When virial 
equilibrium is satisfied, we expect that the specific thermal energy 
in a halo of mass M and of radius R will scale with its potential 
energy, GM/R, while the kinetic energy is proportional to M^^'^. 
Since ctdm expresses the specific thermal energy in dark matter, we 
can express the mean expected velocity dispersion as a function of 
mass as 

CTDM (M, z) = CTDM.IS I -^qis j^^ I ■ (16) 

Here the fit parameters are the slope a of the relation, and the nor- 
malization CTDM.IS at a mass scale of 10^^ h^^ Mq. While the slope 
a just follows from the virial theorem if halos form a roughly self- 
similar family of objects (which they do to good approximation), 
the amplitude ctdm, is of the relationship is a non-trivial outcome 
of numerical simulations and reflects properties of the virialization 



process of the halos as well as their internal structure, fevr ard et al.l 
112008) showed that a single fit is consistent with the numerical data 
of a large set of N-body simulations of the ACDM cosmology, cov- 
ering a substantial dynamic range. 

However it is conceivable that the amplitude of the relation- 
ship will be slightly different in early dark energy cosmologies, as 
a result of the different virial overdensity that is predicted by the 
top hat collapse in these cosmologies. If true, this would then also 
hint at a different normalization of the relationship between total 
Sunyaev-Zeldovich decrement and mass, which would hence di- 
rectly affect observationally accessible probes of the cluster mass 
function at high redshift. 

We here test whether we can find any difference in this re- 
lationship for our different dark energy cosmologies. In Figure [TTl 
we plot the velocity dispersion of halos as a function of mass, in the 
four different cosmologies we simulated. The halos were identified 
using a spherical overdensity definition, where the virial radius r2oo 
was determined as the radius that encloses a fixed multiple of 200 
times the critical density at the redshift z, and M200 being the cor- 
responding enclosed mass. We then determined the best-fit relation 
obtained from our numerical data (red solid lin es). This fit is i n very 
good agreement with the results obtained bv lEvrardet"ai] ( l200a) 
(dotted blue lines), given by ctdm, is = 1082 ± 4.0kms~^ and 
a = 0.3361 ± 0.0026, a value consistent with the viral expectation 
of Of — 1/3. The insets show the residuals about the fit at redshift 
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z = 0. They have a log-normal distribution with a maximum of 6% 
dispersion (for the DECDM model) around the power-law relation. 
The histograms are well fit by a log-normal with zero mean. 

We find that the halos closely follow a single virial relation, 
insensitive of the cosmological parameters, the epoch and also the 
resolution of the simulation. In particular, we do not find any sig- 
nificant differences for the EDE models, instead, the same form of 
the virial relation is preserved across the entire range of mass and 
redshift in the four simulations. The velocity dispersion-mass cor- 
relation hence appears to be global and very robust property of dark 
matter halos which is not affected by different contributions of dark 
energy to the total energy density of the universe. 

This is a reassuring result as it means that also in the case of 
early dark energy, clusters can be studied as a one parameter fam- 
ily and the calibration of dynamical mass estimates from internal 
cluster dynamics does not need to be changed. Differences in the 
normalization should only reflect more or less frequent halo merg- 
ers and inte ractions, which can introduce an additional velocity 
comp onent jEspino-Briones et alJIlOOTi : iFaltenbacher & Mathewsl 
l2007h . 



5.2 Halo concentrations 

As we have seen, for an equal normalization of the present-day lin- 
ear power spectrum, the dark matter halo mass function at 2: = 
does not depend on the nature of dark energy. On one hand this is 
a welcome feature, as it simplifies using the evolution of the mass 
function to probe the expansion history of the universe, but on the 
other hand it disappointingly does not provide for an easy handle 
to tell different evolutions apart based only on the present-day data. 
However, a discrimination between the models may still be made if 
the internal structure of halos is affected by the formation history, 
which would show up for example in their concentration distribu- 
tion. 

Cosmological simulations have consistently shown that the 
spherically averaged mass density profile of equilibrium dark mat- 
ter halos are approximately universal in shape . As a result, we can 
descr i be the halo profiles by the NEW formula jNavarro et aljl993 . 
Il99dll997l) : 

= , (17) 

Pcrit {r/rs){l + r/rsY' 

where pcrit ~ SHq / SivG is the critical density, Sc is the charac- 
teristic density contrast and Vs is the scale radius of the halo. The 
concentration c is defined as the ratio between r2oo and Vs- The 
quantities Sc and c are directly related by 

S _ 200 
' 3 [ln(l + c)-c/(l + c)]- 

The concentration c is the only free parameter in Eqn. dlSb at 
a given halo mass and these two quantities are known to be corre- 
lated. In fact, characteristic halo densities reflect the density of the 
universe at the time the halos formed; the later a halo is assembled, 
the lower is its average concentration. 

We have measured concentrations for our halos in the different 
cosmologies using the same procedur es as applied to the analysis of 
the Millennium simulation jNeto et al. 2007; Gao et al. 2008). Eor 
our measurements, we take into account both relaxed and unrelaxed 
halos. In the second case, the equilibrium state is assessed by means 
of three criteria: (1) the fraction of mass in substructures with cen- 
ters inside the virial radius is small, /sub < 0.1, (2) the normalized 



offset between the center-of-mass of the halo rem and the poten- 
tial minimum Vc is small, s — \rc ~ T'cm|/r2oo < 0.07 and (3) 
the virial ratio is sufficiently close to unity, 2T/\U\ < 1.35. These 
quantities provide a measure for the dynamical state of a halo, and 
considering these three conditions together guarantees in practice 
that a halo is close to an equilibrium configuration, excluding the 
ones with ongoing mergers, or with strong asymmetric configura- 
tions due to massive substructures. 

Eor all relaxed halos selected in this way, we computed a 
spherically averaged density profile by storing the halo mass in 
equally spaced bins in log^Q (r) between the virial radius r2oo and 
logjQ(r/r2oo) ~ —2.5. We used 32 bins for each halo and we 
choose a uniform radial range in units of r2oo for the fitting pro- 
cedure so that all halos are treated equally, regardless of the mass. 
We find that we obtain stable results w hen we use halos w ith more 
than 3000 particles, consistent with the lPower et al. |(|200l) criteria, 
while with fewer particles we notice resolution effects in the con- 
centration measurements, as both the gravitational softening and 
discreteness effects can artificially reduce the concentration. The 
final mass range we explored is hence 10^^ to 10^^ /i~^Mq. 

In Eigure[T2l we show our measured mass-concentration rela- 
tion for the different dark energy models at 2: — 0. The four solid 
lines show the mean concentration as a function of mass. The boxes 
represent the 25 and 75 percentiles of the distribution, while the 
whiskers indicate the 5 and 95 percentiles of the distributions. We 
note that the scatter of the concentration at a given mass is very 
close to a log-normal distribution. It is interesting to remark that 
both the mean and the dispersion decrease with mass. In fact, mas- 
sive halos form in some sense a more homogeneous population, 
because they have collapsed recently and so the formation redshift 
is relatively close to the present epoch. On the other hand, less mas- 
sive halos have a wider distribution of assembly redshifts and the 
structure of individual objects strongly depends on their particular 
accretion histories. Eor them, the assumption that objects we ob- 
serve are just virialized is therefore inappropriate, especially for 
very low mass halos. In Eig. [T2]we take into account only the re- 
laxed halos, but we did an analogue measurement also for the whole 
sample, shown in Eigure[T3lat redshift (top panel) and at z = 1 
(low panel). 

The correlation between mass and concentration approxi- 
mately follows a power law for the relaxed halos of the ACDM 
model. In the literature, the concentrations would be expected to be 
somewhat lower if a complete sample is considered that includes 
disturbed halos. Comparing Eigures [T2I and [T3l we notice that this 
expectation is confirmed, but the difference is not very pronounced, 
only about 5% for the whole mass range. We also note that the 
normalization erg = 0.8 used for our simulat ions slightly lowers 
the amplitude of the zero point of the relation l lMaccio et al.ll200§) 
when compared to the WMAP-3 normalization, as halos tend to 
assemble later with lower as and/or Sim.. 

When we compare our four simulated cosmologies we find 
that, as expected, EDE halos of given mass have always higher 
concentration at a given redshift than models with a cosmological 
constant: they tend to form earlier and so they have a higher char- 
acteristic density. Nevertheless, the differences are not large, they 
deviate by no more than ^ 27% at z = over the entire mass range 
we studied for all halos and ~ 25% for the relaxed one. At higher 
redshift, the differences are only slightly bigger, of order of ~ 28% 
at z = 1 for the whole sample, and ~ 35% for halos in equilibrium 
configuration, suggesting that we anyway need reliable numerical 
calibrations and highly accurate observational data to discriminate 
between the different cosmologies. Interestingly, the average con- 
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Figure 12. Mass-concentration relation for relaxed halos in all our simulations. The boxes represent the 25 and 75 percentiles of the distribution with respect 
to the median value, while the whiskers show the 5 and 95 percentiles. We co mpare our results w ith the theoretical expectations from NFW, ENS, BOl. Also, 
a modified NFW prescriptions with slightly modified parameters as updated bv lGao et alj j2008l) is shown (see Section 5 for details). 



centration is almost independent of mass when we consider z ^ 2, 
as the average concentration of the more massive halos is similar 
at all redshifts jOao et ah 1 12003) and we are then restricted to the 
exponential tail of the mass function. 

The change in concentration normalization relative to the cos- 
mological constant model is well represented by the ratio between 
the linear growth factor of different models at very high redshift, 



Co 



CO,ACDM - 



(19) 



-D+,acdm(oo) ' 

as suggested by iDolag et alj j2004l) . In Table |2] we compare the 
ratio between the concentration at z = both for the relaxed ha- 
los (second column) and for the whole sample (third column) with 
the ratio between the asymptotic growth factor for the same cos- 
mologies (forth column). The order of magnitude of the two effects 
is comparable, although the match is not perfect. Here the ratios 
are computed for i\f ~ 4 x 10^^ Ii^^Mq, where we have a large 
number density of halos. 

It is interesting to compare the concentrations we measure 
with the various theoretical predictions that have been made for 
this quantity. We investigate three popular descriptions for the con- 
centration: the class ic Navarro, Frenk & White model (hereafter 
N FW), the model o f ISullock et alj j200lL hereafter BOl), and that 
of lEke et"aL I ( I2OOII, hereafter ENS). Finally, we also plot the new 
modified version of the original Navarro Frenk and White model, 
as recently proposed bv lGao et alj ^20081) . Both the BOl model and 



Model 


CO 




13+ (inf) 


CO, ACDM 


f^Q.ALL.ACDM 


-D+,ACDM(int) 


ACDM 


1.000 


1.000 


1.000 


DECDM 


1.256 


1.275 


1.228 


EDEl 


1.218 


1.232 


1.229 


EDE2 


1.255 


1.273 


1.252 



Table 2. Concentration and asymptotic growth factor in the four different 
cosmologies studied here. The ratio between the concentration parameters 
at redshift z = refers to a mass of M ~ 4 X 10^^ h~^MQ that corre- 
sponds to the mass range that contains the majority of our halos. For each 
model (first column) we give the cq parameter relative to ACDM, taking 
into account the relaxed halos (second column) or the whole sample (third 
column). Finally, in the last column we show the linear growth factor at 
infinity relative to the ACDM cosmology. 



the standard NFW have two free parameters that have been tuned to 
reproduce simulation results. In the original NFW prescription, the 
definition of the formation time of a halo is taken to be the redshift 
at which half of its mass is first contained in a single progenitor: 
F = 0.5. The second parameter is the proportionality constant, 
C = 3000, that relates the halo density scal e to the mean co smic 
density at the collapse redshift Zcoii- Recentlv. lGao et al no- 
ticed that the evolution of the mass-concentration relation with red- 
shift can be approximated much better by setting F = 0.1. The 
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Figure 13. Mass-concentration relation for all halos in our simulations. The 
top panels refers to redshift z = 0, while the bottom panel shows the results 
at 2 = 1. The boxes represent the 25 and 75 percentiles of the distribution 
with respect to the median value, while the whiskers show the 5 and 95 
percentiles. We compare our results with the theoretical expectations from 
NFW, ENS, BOl. Also, a mo dified NFW pres cription with slighdy modi- 
fied parameters as updated bv lGao et alj j2008l) is shown (see Section 5 for 
details). The concentration is 5% lower with respect to the relaxed sample 
at 2 = for the ACDM model. 



BOl model adopts as collapse redshift the epoch at which the typ- 
ical collapsing mass fulfills A/, (ac) = F Mvm, with F = 0.01. 
They further assume that the concentration is a factor K = 3.4 
times the ratio between the scale factor at the time the halo is iden- 
tified and the collapse time. For K and F we use the values that 
are indicated as the best parameters by iMaccio et alj j2007h . Fi- 
nally, we compute the ENS prescriptions considering the effective 
amplitude of the power spectrum at the scale of the cluster mass. 
Tiiis quantity, rescaled for the linear growth factor of the simulated 
cosmology, has to be constant. In this cas e, onl y one parameter , 
Ca = 28, is needed. iBuUock et al] ('2001') and E ke et alj j200lh 
refer to the virial radius as the one including an overdensity given 
by the generalized top hat collapse model. We have appropriately 
adapted these models such that the concentration of a halo is de- 
fined instead relative to radius r2oo, as in the NFW model. 

Aside from BOl, all three other model predictions yield con- 



centrations that agree reasonably well with the measured values 
at 2: = 0. The BOl model underpredicts the relation at high 
masses, where it gives has a sharp decline of the relation for 
M 10^^ /i~^M0 which is not seen in the simulations. In contrast, 
the NFW model is in reasonable agreement with the data at z = 
for both halo samples. However, at z = 1 the evolution predicted by 
the NFW model is less than what we find numer ically, even when 
we consider the revised formulation proposed bv lGao et alj j2008l) 
(indicated as NFW modified). The NFW model with the new fit- 
ting parameters yields a reasonable fit at the high mass end, but 
performs a bit worse than the original formulation at z = 0, spe- 
cially at low masses. Unfortunately, for the NFW model the nor- 
malization is model dependent, so we cannot really capture all the 
effects due to different cosmological parameters we use. Finally, 
the dashed black line in each plot shows the ENS model. This pre- 
scription gives the best match with our results and has been able 
to reproduce the slope of the concentration-mass relation even at 
higher redshift. 

At a fixed mass, halos in the EDE cosmology are significantly 
less concentrated than their counterparts in the ACDM cosmology. 
It is interesting to notice that the ENS model reproduces these dif- 
ferences quite well, without modifications of the original prescrip- 
tion. In Figure [141 we plot for each simulation the corresponding 
theoretical expectation (dashed lines) for the sample of relaxed ha- 
los at z = 0. For a low density universe the scaling of the linear 
growth factor with redshift leads to a greater difference between 
the models. Dark halo concentrations depend both on the redshift 
evolution of 5c and the amplitude of the power spectrum on mass 
scales characteristic for the halo. 

These results for the concentration are particularly impor- 
tant since they demonstrate that quintessence cosmologies with the 
same equation-of-state at present, but different redshift evolution, 
can produce measurable differences in the properties of the non- 
linear central regions of cluster-sized halos. However, the prospects 
to observationally exploit these concentration differences to distin- 
guish different dark energy cosmologies are sobering. For one, the 
systematic differences we measure for the concentrations are quite 
small compared to the statistical errors for the mean concentration, 
while at the same time the theoretical algorithms for predicting the 
halo concentration perform quite differently already for the ACDM 
cosmology. Furthermore, directly measuring halo concentrations in 
observations is not readily possible as it requires an accurate knowl- 
edge of the virial radius of a halo, a parameter which is poorly con- 
strained from observations. It therefore remains to be seen whether 
the effects of dark energy on the non-linear structure of dark halos 
can be turned into a powerful tool to learn about the nature of dark 
energy. 



6 COUNTING HALOS BY VELOCITY DISPERSION 

As we have seen, the different evolution of the halo mass function 
is in principle a very sensitive probe of the expansion history of 
the universe, especially when the massive end of the mass function 
is probed. Obtaining absolute mass estimates from observations is 
however problematic, and fraught with systematic biases and un- 
certainties. It is therefore important to look for new ways to count 
halos which are more readily accessible by observations. 

One such approach lies in using the motion of galaxies in 
groups or clusters of galaxies to measure the line-of sight veloc- 
ity dispersion, which in turn can be cast into an estimate of the total 
virial mass of the host halo. This relies on the assumption that the 
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Figure 14. Mass-concentration relation for relaxed halos today. Here we 
show the agreement between simulation results (symbols) and theoretical 
predictions from ENS (dashed lines), both for the ACDM and EDE cos- 
mologies. To this end we solve Eqn. (13) and (16) of Eke et al. (2001). The 
differences between the four cosmologies are due mostly to the differences 
in the growth factor evolution and consequently in the amplitude of the 
power spectrum. The ENS formula works quite well also for EDE models 
without modifications of the original prescription. 




Figure 15. The velocity function n(cr) as a function of halo mass, for all 
satellites inside r200 ■ The shaded area indicates the differences between a 
ACDM model with eg = 0.8 and the same model with erg = 0.9. It is 
interesting to remark that this EDE models could justify a higher normal- 
ization cosmology. 



dynamics of the cluster or group galaxies is tracing out the dark 
matter halo potential. 

Cluster and group galaxies can be ide ntified with dark mat- 
ter sub-structures in N-body simulations dSpringel etai] 1200 ll : 
IVale & Ostrike3l2004l) . Employing the bulk velocities of sub-halos 
as a simulation proxy for real galaxy velocities, we can hence build 
a velocity profile for any isolated halo, and estimate a line-of-sight 
velocity dispersion, similarly as it is done for observed group cat- 
alogues of galaxies. This allows then to directly count halos (i.e. 
galaxy groups) as a function of line-of-sight velocity dispersion, 
bypassing the problematic point of assigning halo mass estimates. 

In Figure [15] we show our estimated cumulative velocity 



dispersion function for our four different cosmologies at redshift 
z — 1.5. This graph can be interpreted as being a different rep- 
resentation of the halo mass function, except that it is in principle 
directly accessible by observations. For this measurement, we have 
derived the information on the velocities from the SUBFIND algo- 
rithm directly implemented in GADGET-3, which can find subhalos 
embedded in dark matter halos. 

An important aspect of this statistic is that it does not rely 
on the often ambiguous definition of a group mass. Instead, it can 
be directly measured and is more readily access ed by observa- 
tions. In fact, studies based on th e DEEP2 survey dLin et alj|2004 
iDavis et al. I l2005l : IConrov et alj [2007) indicate that, if combined 
with both the velocity dispersion distribution of clusters from the 
Sloan Digital Sky Survey and independent measurements of as, 
they will be able to constrain w to within approximately 1% accu- 
racy. This method is almost independent of cosmological parame- 
ters, with the exception of ag, since a change in normalization can 
shift the space density of halos as a function of mass by a similar 
amount as done by the EDE models. This is illustrated in Figure [Tsl 
by the shaded area, which represents the change of the velocity dis- 
persion function when ag is increased from 0.8 (green line) to 0.9 
(upper limit of the shaded area). The velocity distribution function 
of the EDE models then approaches the one that we would measure 
for a ACDM model with higher erg. 

These kind of studies have strong motivations both from the 
observational and theoretical point of view: there is little scatter 
between host galaxy luminosity and dark matter halo virial mass 
and the velocity difference distribution of satellites and interlop- 
ers can be modeled as a Gaus sian and a constan t, respectively 
jConrov et al. 2005; Faltenbache r & Diemandll2006l) . 

Figure [16] (left panel) shows the cumulative number of groups 
with velocity dispersion above a given value, as a function in red- 
shift for the different models. We decided to count halos above a 
velocity dispersion of 300kms~^, where accurate measurements 
can be expected also from observations. Note that there is already 
a very large difference between ACDM and EDE at redshift z = 1. 
We find that there is almost no evolution in the cluster number in 
the dark energy models, while ACDM drops by a factor of nearly 
10 up to redshift z — 3. What is especially important here is the 
relative difference between the number counts of the two different 
cosmologies. The fact that we do not need to introduce the mass in 
this comparison give us the advantage of having no error derived 
from the particular measurement procedure adopted for the mass. 
At a fixed velocity dispersion, we can directly probe the growth 
of the structure at each redshift, which depends on the equation of 
state parameter w. The slower evolution of the cluster population in 
EDE models is exactly what is expected to be observed also from 
Sunyaev-Zeldovich studies of large samples of clusters of galaxies. 
Combined with probes of the cluster internal velocity dispersion 
we can hence hope to be able to derive stringent cosmological con- 
straints. 

We also remark that the relative difference between the num- 
ber of objects within these four simulations seems to be a quite 
robust statistic which is invariant with respect to details of the mea- 
surement procedure. For example, in Figure [16] (right panel), we 
change the number of considered subhalos in the halos to be a min- 
imum of 3, 4, and 5, but the velocity dispersion function relative 
to the ACDM cosmology is essentially unchanged. In practice, the 
number of observable satellites per host halo suffers from limita- 
tions imposed by the magnitude limit of the survey. Our results 
suggest that the measured velocity dispersion should however be 
relatively insensitive to this selection effect. 
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Figure 16. Left panel: Comparison of the redshift evolution of the velocity dispersion function for all four cosmologies we simulated (ACDM, DECDM, 
EDEl, and EDE2). Here the cumulative count of groups with velocity dispersion above cr = 300 kms~^ was used to measure the amphtude of the velocity 
dispersion function. Right panel: Differences in the number count when only halos with more than 3 (solid line), 4 (dotted line) or 5 (dashed line) substructures 
are selected. 



Finally, we have also studied a few properties of the largest 
substructures in halos to see whether there is a difference in EDE 
cosmologies. In Figure [771 the small diamonds indicate the values 
of the ratio between Mi (the mass of the most massive subhalo) 
and M200 (the mass within a sphere of density 200 times the critical 
value at redshift 0) for the first 200 most massive halos at redshift 
2 = 3. The filled circles represent the median of the distribution, 
computed in bins of 50 halos each, while the error bars mark the 
20-th and 80-th percentiles of the distribution. There is almost no 
dependence on parent halo mass, but we can notice a small, but 
systematic tendency for the ACDM subhalos to be slightly more 
massive. The dependence is quite weak, yet this behaviour is clear 
even if the mass of the progenitor A/200 tends to be lower on av- 
erage at this high redshift. This is symptomatic of the fact that the 
ACDM substructures are formed at lower redshift with respect to 
what happens in the EDE models. This is also consistent with ex- 
pectations based on the obse rved dependence of s ubstructure mass 
fraction on halo mass (e.g. iDe Lucia eta"i]|2004l) . Once accreted 
onto a massive halo, subhalos suffer significant stripping, an effect 
that is more important for substructures accreted at higher redshift, 
making the subhalos in the EDE models less massive on average. 
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Figure 17. Ratio of the mass of the two most massive substructures with re- 
spect to the mass of the parent halo. The small diamonds refer to individual 
halos, while the filled ckcles are the median values. The en'or bars mark the 
20th and 80th percentiles of the distributions. 



7 CONCLUSIONS 

In this study we have analyzed non-linear structure formation in 
a particular class of dark energy cosmologies, so called early dark 
energy models where the contribution of dark energy to the total en- 
ergy density of the universe does not vanish even at high redshift, 
unlike in models with a cosmological constant and many other sim- 
ple quintessence scenarios. Our particular interest has been to test 
whether analytic predictions for the halo mass function still reliably 
work in such cosmologies. As the evolution of the mass function is 
one of the most sensitive probes available for dark energy, this is 
of crucial importance for the interpretation of future large galaxy 
cluster surveys at high redshift. The mass function of EDE mod- 
els is also especially interesting because analytic theory based on 
extensions of the spherical collapse mod el predicts that the mas s 
function should be significantly modified jSartelmann et alj2006h . 



and in particular be characterized by a different value of the lin- 
ear overdensity 5c for collapse, as well as a slightly modified virial 
overdensity. 

We have carried out a set of high-resolution N-body simu- 
lations of two EDE models, and compared them with a standard 
ACDM cosmology, and a model with a constant equation of state 
equal io w = —0.6. Interestingly, we find that the universality of 
the standard Sheth & Tormen formalism for estimating the halo 
mass function also extends to the EDE models, at least at the 
^ 15% accuracy level that is reached also for the ordinary ACDM 
model. This means that we have found good agreement of the stan- 
dard ST estimate of the abundance of DM halos with our numer- 
ical results for the EDE cosmologies, without modification of the 
assumed virial overdensity and the linear density contrast thresh- 
old. This disagrees with the theoretical suggestions based on the 
generalized top-hat collapse. In fact, if we instead use the latter as 
theoretical prediction of the halo mass function, the deviations be- 
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tween the prediction and thie numerical results become significantly 
larger. We hence conclude that the constant standard value for the 
linearly extrapolated density contrast can be used also for an analy- 
sis of early dark energ y cosmologies. Very recently, similar results 
were also obtained bv lFrancis et alj ( 1200 Sh . who studied the same 
problem in cosmological simulations with somewhat smaller mass 
resolution. 

This results on the mass function appear to hold over the 
whole redshift range we studied, from z = to 2 = 3. Since our 
simulations were normalized to the same as today, their mass func- 
tions and power spectra agree very well today, but towards higher 
redshift there are significant differences, as expected due to the dif- 
ferent histories of the linear growth factor in the different cosmolo- 
gies. In general, structure in the EDE cosmologies has to form sig- 
nificantly earlier than in ACDM to arrive at the same abundance 
today. For example, already by redshift z = the abundance of 
galaxy clusters of mass A/ = 5 X 10^^ /i-^Mq is higher in EDEl 
by a factor of '--^ 1.7 relative to ACDM. 

The earlier formation of halos in EDE models is also directly 
reflected in the concentration of halos. While for a given as we find 
the same abundance of DM halos, the different formation histories 
are still reflected in a subtle modification of the internal structure of 
halos, making EDE concentrations for all halo masses and redshifts 
considered slightly higher. The difference is however quite small, 
but it would, for example, lead to a higher rate of dark matter anni- 
hilation in halos. 

Another relationship that appears to accurately hold equally 
well in ACDM as in generalized dark energy cosmologies is the 
viri al scaling between mass and dark matter velocity dispersion 
that lEvrardetai] j2008l) has found. In fact, we find that their nor- 
malization of this relation is accurately reproduced by all of our 
simulations within the measurement uncertainties, independent of 
cosmology. This also suggests that possible differences in the virial 
overdensity of EDE halos must be very small, and that presumably 
the relationship between total Sunyaev-Zeldovich decrement and 
halo mass is unmodified as well. 

We show that counting the number of halos as a function of 
the line-of-sight velocity dispersion (of subhalos or galaxies), both 
in simulations and observations, can probe the growth of struc- 
tures with redshift, and so put powerful constraints on the equa- 
tion of state parameters. This goal can be achieved by just identi- 
fying and counting groups in galaxy survey data such as DEEP2, 
and by comparing them with high-resolution N-body simulations. 
Precision measurements with this technique will still require accu- 
rate calibrations to deal with complications such as a possible ve- 
locity bias or selec tion effects in observational surveys. However, 
iDavis et al. I j2005l) suggest that the DEEP2 survey alone has the 
power to constrain w to an accuracy of 20% using velocity disper- 
sion data, which illustrates the promise of this technique. In combi- 
nation with other independent data, such as X-ray temperature and 
SZ decrement data, the constraints could be improved to an accu- 
racy of 5%, without the need to invoke a model for the ambiguous 
total mass of a halo. 

Distinguishing a time-varying dark energy component from 
the cosmological constant is a major quest of the present theoretical 
and observational astronomy. One approach is to rely on classical 
cosmological tests of the Hubble diagram, e.g. by pushing the su- 
pernova type la observations to much higher redshift. Another quite 
direct geometrical probe is the observation of baryonic acoustic os- 
cillations in the matter distribution at different redshifts. Finally, 
the linear and non-linear evolution of cosmic structures provides 
another opportunity to constrain dark energy. In this work we have 



used numerical N-body simulations to examine the difference in 
structure growth in early dark energy cosmologies. We have seen 
that such simulations are essential to test the predictions of more 
simplified analytic models, and to calibrate observational tests that 
try to constrain the properties of dark energy with the abundance 
and internal structure of dark matter halos. Our results show clearly 
that the effects due to dynamical dark energy tend to be quite sub- 
tle, and can only be cleanly distinguished from ordinary ACDM in 
high accuracy simulations. This poses new challenges to improve 
the precision of future generations of simulations, and at the same 
time emphasizes the immense observational task to arrive at suffi- 
ciently precise data at high redshift to constrain the dark side of the 
universe with the required accuracy. 
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